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1. Introduction 



The study of B-physics plays a fundamental role within flavour physics both in accurately 
testing the Standard Model and in the search of New Physics effects. To this aim it is crucial to have 
theoretical uncertainties under control, in particular those of the hadronic parameters computed on 
the lattice. 

With the available computer power it is not possible to simulate quark masses in the range 
of the physical b mass keeping, at the same time, finite volume and discretisation effects under 
control. In order to circumvent these problems, many different methods have been proposed so far 
(see ref. [p for an up to date collection of results). 

The approach that we have adopted and that we discuss below consists in using lattice QCD 
data with the heavy quark mass ranging from the charm region up to ~ 4/5 of the physical b 
quark mass, together with the information coming from the static limit point. In order to deal with 
the simulated light quark mass and finite lattice spacing, a careful extrapolation to the chiral and 
continuum limits has been performed. An alternative method, based on the introduction of suitable 
ratios having an exactly known static limit, has been recently proposed and investigated by our 
Collaboration (ETMC) [g]. 

In section || we describe the computation of the decay constants in the static limit; in section [| 
we present the interpolation between the charm and infinite-mass sectors and compare the results 
with those obtained in ref. [g]. 



2. Heavy-light decay constant in the static limit of HQET 

We have combined a light doublet of twisted-mass fermions (\j/ T = (u,d)) defined at maximal 
twist with a static quark described by the HYP2 action [Bl] to improve the signal-to-noise ratio [0J : 



= a 4 £ Yh (x)V* y h (x) , V* y h (x) = - y h (x) - U^ YY2 (x-aO) y h (x-aQ) 



(2.1) 



In order to extract the decay constant using maximally twisted lattice QCD, we need to evaluate the 
matrix element of the static-light local current. At maximal twist the pseudoscalar current (£P stat ) R 
in the physical basis, in terms of the twisted basis used in the numerical simulations (light quark 
fields x T = (Xu,Xd)), is given by 

(^ sta % (x) = (%(x) Y5 u(x)) R = ±= (Z$«P(x) + iZ$*S(x)) (2.2) 

where P = YhJsXu an d S = YhXu are th e pseudoscalar and scalar densities which renormalise with 
the Zp tat and Z| tat appropriate to the static-light framework. 

We define c\ = i {0\tyh% u \B) and C5 = (0|i//h75Z«|#) where \B) is the lattice ground state. At 
maximal twist, the amplitude we need to compute is <I> = /bVMb = (2| tat ci +Z s p at cs). The (bare) 
matrix elements c\ and C5 have been measured from an analysis following the static HQET spec- 
trum study with twisted-mass fermions [Qj. The ETMC ensembles 61,2,3,4 and C1.2 [§, g] have so 
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far been considered (i.e. two lattice spacings). Here we concentrate on the lightest heavy-light 
meson state, the pseudoscalar meson which we call here the B meson (or B s with a strange valence 
quark). We take the value of m q for the strange quark from the ETMC studies of the strange-light 
mesons [||, ^] which used the same gauge configurations as used here, namely am s = 0.022 at 
j3 = 3.9 and 0.017 at j8 = 4.05. We measure the correlation of operators at source and sink with a 
large choice of operators: local and smeared; parity conserving and non-conserving. We then make 
a simultaneous fit to a sub-matrix (typically 6 x 6 ) in a given Euclidean time t interval. We chose 
this ^-interval to have similar physical extent at different lattice spacings. We find that the non-local 
operators have weaker coupling to excited states, as expected. Such non-local operators can give a 
good determination of the energy levels but to extract the required matrix element (related to fp ) 
we need to include local operators in the fit. At j3 = 3.9 we use a 4 state fit with t/a range 4 — 10 
but with the correlations that have local operators (at sink and/or source) we restrict to t/a range to 

6 — 10. This choice gives acceptable values of % 2 using correlated fits. We then make uncorrected 
fits to determine the required energies and matrix elements with statistical errors determined by 
bootstrap. At jS = 4.05 the appropriate t/a range is found to be 5 — 12 for smeared correlators and 

7 — 12 for local ones. We have checked by making many different fits that the fit parameters are 
stable, within the statistical error assigned. For the correlations of B s mesons, we make similar 
fits but find that the minimum t/a value has to be increased by 1 unit to preserve an acceptable 
(correlated) % 2 . 

Then one computes Z), tat and Z| tat in order to get the matrix element renormalised in HQET 
at a specific scale pt. We have chosen to renormalise it in the MS scheme at /I = I /a and for this 
preliminary account of our work the renormalisation is done perturbatively at 1 loop order. MS is a 
continuum-like scheme defined within dimensional regularisation, while the regulator of our bare 
quantities is the inverse lattice spacing. So one needs a matching between both regularisations. It 
can be written as 



(0(p,ii)) DRMS = [l-^(-%m«V + C°) (0( P ,a)) 



lat 



a s 

Z (all)(0(p,a)) M , (2.3) 



where the renormalisation scheme and scale of the coupling constant a s is not specified at this level 
of perturbation theory. Expressions of C p ^ are complicated and not illuminating, essentially due 



to the HYP-smeared static action and the improved part of the gluon propagator [10]. Thus we 
have simply collected the numerical values of Zf, tat and Z| tat in Table [j] for a boosted coupling gp = 
Sq/ {Up} (where = 6/j3 and (Up} is the average plaquette value). It turns out that the systematic 
error introduced by a poor determination of the ratio z r = Zf at /Z p tat is minimal, especially on the 



p 


-vstat 


ystat 


3.9 
4.05 


0.849 
0.859 


0.933 
0.938 



Table 1: First order perturbation theory renormalisation factors of the pseudoscalar and scalar static-light 
dimension 3 operators in the MS scheme at the scale \l = \/a. 
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Figure 1: Left plot: unrenormalised heavy-light decay constant combination r ' <f>/Z(4.05) (with Z = 
(Z|, tat + Z| tat ) /2) versus the squared mass of the pion built of the light sea quarks. The circles represent 
the B meson case, where the valence light quark is equal to the sea quark. The squares represent the B s 
meson case, where the light valence quark is the strange quark. The data at j3 =3.9 (red symbols) have been 
multiplied by the appropriate factor to match the same scale for the data at j3 = 4.05. The curves represent 
the NLO HMChPT theory expressions. Right plot: the ratio versus the squared mass of the pion built of 
the light sea quarks. The curve represents the NLO heavy quark chiral perturbation theory. 



ratio of the B and B s decay constants. We thus present in fig. |] the bare matrix element, which 
depends on the ratio z r only. 

Once the matrix element <J> M5 (/i = \/a) has been renormalised in the MS scheme at the scale 
|l = l/aa NLO running of perturbation theory [ LI ] has been applied to evolve it to a scale jU = 
M^ xp . This is what is needed to perform a fit together with the relativistic data matched to HQET 
at the same scale (see next section). 

The extrapolation of <J># down to the physical pion has been performed with Heavy Meson 



Chiral Perturbation Theory (HMChPT) at NLO by using the formula []L2J-|14|] 



<I>o 



3(1 + 3| 



a2\ 



Ml 



(4nfy 



■log 



{Anff 



+ a l Mf l . 



(2.4) 



where Mu denotes the simulated pion masses, / stands for the light decay constant in the chiral 
limit, while Q>q( s ) and a,u s ) are free fit parameters. The g 2 coupling has been fixed to 0.2 [JTsJ, p| ], 
and we have checked that a change of 50% in the value of g 2 results in a shift in <I>b which is well 
below the statistical error. The chiral extrapolation of <J>g, <t>s s and the ratio QbJ'&b is shown in 
fig. |. This fi gure also illustrates that we find consistent results at our two available lattice spacings 
within the relatively large errors. We do not have enough data to include explicit discretisation 
error terms in the fit formula. However it seems that cut-off effects are quite small. This is more 
evident for the ratio which is consistent with having no cutoff effects (see right plot of fig. [T]). 
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P=4.05 (a = 0.065 fm), V/a 4 = 32 3 x 64 





■ 


3*1^=0.2000, a*(l = 


=0.0030 




• 


a*|l h =0.2000, a*)l s 


=0.0220 


1,3 




a*|l h =0.4190, a*|i = 


=0.0030 




♦ 


a*|l h =0.4190, a*n_ 


=0.0220 




Figure 2: Effective masses at j3 = 4.05 for two heavy-light (hi) and two heavy-strange (hs) quark mass 
combinations. The two heavy quark masses correspond approximately to the physical charm quark mass 
and to ~ 2/3 of the value of the physical b quark mass. 



3. Relativistic results and interpolation to the physical b quark mass 

We perform an interpolation of the heavy-light (hi) decay constants from the charm region up 
to the bottom mass, by including data in the static limit calculated in the HQET as explained in 
the previous section. The lattice QCD data used in this analysis are at four values of the lattice 
spacing a 0.100, 0.085, 0.065, 0.050fm (corresponding to j8 = 3.8,3.9,4.05,4.2), that is we have 
used the configuration ensembles denoted in []6[ ^] as A2.3, Bi,2,3,4,6,7> Ci.2,3 and D2, respectively. 
We have simulated for each ensemble 16 heavy quark masses in the range m^ xys < mj, < Q.8m? hys . 
Quark propagators with different valence masses are obtained using the so called multiple mass 
solver method [|l7|]. In fig. ||| we show for illustrative purpose the effective masses at j8 = 4.05 and 
for few quark mass combinations. 

The analysis is performed by studying the dependence of the decay constants, more precisely 
of the quantity <J>/^ = fhq\fMhq, as a function of the meson masses, as in our recent analysis of the 
fo and fo s decay constants [§J. 

In order to make use of the HQET scaling low we introduce for each simulated hq meson mass 



Mh q the HQET quantity that is finite in the static limit [11] 



0>, 



hq 



a™S(M hq ) \ 
MS( M - P) J 

m (M hq ) 



-»/(2ft) 



a 



439 



28?r 2 \ a MS (Mhq) - a MS (M| xp ) 



1089 297 



4k 



1 + 



8a r 



An 



(®hq)QCD , 



(3.1) 



which has been obtained through the NLO matching from QCD to HQET and evolving at NLO to 
the renormalisation scale given by the experimental value of the B meson mass. For <J>/, ? (q = l,s) 
we first study the dependence on the light/strange quark mass at fixed heavy mass through the 
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following functional forms 

= A'(a,m h ) ■ (1 +B' ■ M?, + C'(a) ■ M 2 S ) . (3.2) 



We note that the fit forms above follow from the HMChPT formulae [12-14], which we have al 



ready used in the static sector (see eq. (2.4)). A dependence of the coefficients A,A',C' on the lattice 
spacings is allowed, in order to account for discretisation effects. The extrapolation/interpolation 
to the physical light/strange quark mass is performed by replacing in eq. ( |3.2| ) M 2 , = (M^ xp ) 2 , 
M 2 S = 2 (M e ^ p ) 2 — (M^ p ) 2 . This first step provides the values of the decay constants at the physical 
light/strange quark mass for every simulated lattice spacing and heavy quark mass, or equivalently 
the quantities <& hqP h ys . 

The second step consists in studying the dependence of <& hqP hy S , included the available static 
points, on the heavy quark mass and on the lattice spacing, in order to interpolate to the b quark 
mass and to extrapolate to the continuum limit. Several functional forms with different £?(a 2 ) and 
&{a ) discretisation terms have been tried, which can be written in a compact way as 

= Y, p nk a2 " M t n q ~ k ' (« = 0, 1,2; * = 0,1,2) , (3.3) 

n,k 

where M/ lq is a reference meson mass with the same simulated heavy quark mass as in the fitted 
quantity <1> and the light quark mass is fixed to a similar value for all data. We have performed 
correlated fits by assuming the static results uncorrected with the relativistic data. 



The results for the decay constants f B and fg s are finally obtained by replacing in eq. ( p.3| ) 
Mf, q = = M e B xp , setting the lattice spacing equal to zero and performing the matching from 
HQET back to QCD at NLO. 

The dependence of the decay constants on the hq meson mass is shown in fig. [| where, for 
illustrative purpose, we also show curves corresponding to one of the various fits. The discretisation 
terms included in the shown fits are of &{a 2 M hq ), Gi^M^ and 6{a A M\ q ) for both ^ MP h ys and 
^hspfy* ■ W e observe that with our data it is not possible to determine the coefficients of more than 
three discretisation terms for each fit and that, in some cases, only two out of three parameters turn 
out to be different from zero. About twenty of these fits have a chi square per degree of freedom 
of order one or smaller and are considered in deriving our final result for fs and fg s . The spread 
among these fits is included in the systematic uncertainty. 

Our preliminary results for f B , f Bs and the ratio read 1 

f B = 191(6)(12)(3)MeV = 191(14)MeV, 
f B , = 243(6)(12)(3)MeV = 243(14) MeV, 
f B Jf B = 1.27(3)(4) = 1.27(5), (3.4) 

where: i) the first error is of statistical plus fitting origin, ii) the second error, estimated through 
the spread of the results obtained with functional forms containing different discretisation terms, 



1 The results given in the present proceedings are based on a larger statistical sample w.r.t to the values presented at 
the Conference and cited in ref. [Q]. 
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Figure 3: Interpolation to the b quark mass and continuum extrapolation of ® M phys (left) and ^ hsP h ys (right). 



represents the residual uncertainty due to the continuum limit and to the b mass interpolation, iii) 
the third error takes into account the effect of the systematic uncertainty on the static point. 



We conclude by comparing the results in eq. ( |3.4| ) with those obtained in ref. [|2p using suitable 
ratios having an exactly known static limit. The latter values read 



f B = 194(16)MeV, 
f Bs = 235(ll)MeV, 



(3.5) 



where the uncertainty is the sum in quadrature of the statistical and systematic errors. The two sets 
of results are in very good agreement, thus providing further confidence on their robustness. We 



note that the results in eq. (3.5) are obtained from a subset of the data analysed in the present study. 
The inclusion of the full set of data is in program for a forthcoming publication. 
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